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I. INTRODUCTION 


Ultra-fast physical processes captured the attention of the scientihc community 
because they arise in different physical situations ranging from chemistry to 
condensed matter Q, to bio physics j^. From the theoretical point of view these 
processes are described by open systems, i.e. systems interacting with the sur¬ 
rounding environment. Most commonly, for an effective description open quantum 
systems are approximated by Markovian dynamics, which require a large separation 
of time-scales between the system and the environment J|. Ultra-fast processes are 
those for which the time-scale of the relevant system is about of the same order 
as that of the bath into which it is immersed. Therefore an effective characteriza¬ 
tion of such processes - if existing at all - must rely on other modes of description, 
which most likely need to be non-Markovian. This justihes the growing interest in 
non-Markovian open quantum system dynamics. It is important to stress that non- 
Markovian descriptions can only be effective on short time-scales: this feature will 
play a crucial role in choosing the best strategy to tackle non-Markovian dynamics. 

The model most widely used to describe open quantum systems is the “indepen¬ 


dent oscillators” (10) model 


7j. In this model the system is bilinearly coupled to 


a bath of independent harmonic oscillators. This model has been thoroughly stud¬ 


ied, and it proved fundamental 

nn 

quantum Brownian motion p, [8 


or the description of Markovian and non-Markovian 
9j. Considering the Heisenberg equations of motion 
of this model, one can derive a Generalized Langevin Equation (GLE), that gives a 
phenomenological description of how the environment affects the system |5l. Iinl-ll4|. 
The structure of the GLE is perfect for deriving the Markovian limit, which in the 
classical regime recovers the Langevin equation for Brownian motion. In short, the 
GLE is suitable to analyze thermal effects, diffusive effects, and the fluctuation- 
dissipation relations due to the environment. However, as we shall argue the GLE is 
not suited in situations when the timescale separation is small and non-Markovian 
features are dominant. Indeed, the picture given by the 10 model and the GLE 
only tells us that the interaction with the environment gives rise to some non- 
Markovian/memory effects, but it is not clear how one can capture them in an 
effective way. 
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In order to tackle this problem, we will consider another representation which suits 
better the analysis of short time dynamics. As we shall show the 10 model indeed 
is not suitable for this scope, because the influence of the interaction between the 
system and the environment is not “time’-ordered as all the environmental oscillators 
act “at the same time” on the system. A better representation from this point 
of view is a chain model, i.e. a system interacting with a chain of first neighbor 
interacting harmonic oscillators |15l. Il6| . This model allows to order the influence 
of the interaction in time, and one can have a clearer physical picture of what 
is going on; the hrst oscillator of the chain will first affect the system, and only 
at a later time the influence of the second oscillator of the chain will reach the 
system, and so on and so forth. As we will show, such a model allows to express 
the system dynamics in terms of the microscopic motion of the bath constituents 
[cf. Eq. fl^Uj) ]. As 10 models, due to their structural closeness to GLE are mostly 
considered in applications, we ask the question; How can we achieve a chain structure 
from an underlying model of independent oscillators? The neatest way to perform 
the transformation between these two descriptions is by solving an Inverse Eigenvalue 
Problem (lEP). These problems are well known in the held of vibration theory, and 


the literature is vast (see e.g. {^, 18| and references therein). 


Our analysis applies both to Classical Mechanics (if one consider phase-space 
variables) and to Quantum Mechanics (if one considers Heisenberg equations of 
motion). The paper is organized as follows; in Sec. 2 we introduce the two models 
for the environment, and we solve the lEP that univocally determines the parameters 
of the chain. In Section 3 we solve the equations of motion of the chain model in 
terms of the environmental modes, obtaining the exact dynamics of the system. We 
eventually compare our results with the phenomenological GLE. 


II. INDEPENDENT OSCILLATORS AND CHAIN MODELS 

We consider an open quantum system made of a particle bilinearly interacting 
with an environment of N independent harmonic oscillators. We assume the environ¬ 
ment to be arbitrarily big but hnite. The Hamiltonian describing this “independent 
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oscillators” (10) model reads: 


H,o = 


p 


+ V{x) 


N 

X ^ CkQk 
k=l 


N 


5 ^ ^ (Pfc + ^kQk) , 


( 1 ) 


k=l 


where x,p are the position and momentum operators of the relevant system, 
are positive constants, and qk,Pk are position and momentum operators of the en¬ 
vironmental oscillators with proper frequency Uk- V{x) is a generic renormalized 
potential, which includes the term x^'Yl,k^\/‘^^k guarantees boundedness and 
translation invariance of the model. Note that this implies that also an initially free 
particle will always acquire an oscillatory behavior due to the interaction with the 
environment. Introducing the vector of the position operators of the environmental 
oscillators = (gi,... g^r), one can write their Heisenberg equations of motion as 
follows: 




q(f) = -n; ■ q(f), 


( 2 ) 


where u: is the diagonal matrix of the oscillators frequencies: to = diag(a;f,.. .cj^), 
with uji < 0 J 2 < . ■ -cjn- Substituting the system ([2]) in the Heisenberg equation for 
X, one obtains the following GLE p, l6|: 


x{t) + / dsp{t — s)x{s) + V'{x) = g{t), 


( 3 ) 


where the prime denotes differentiation with respect to x. The kernel p{t — s) is 
called friction kernel, since in the Markov limit the integral term becomes the friction 
constant of the Langevin equation. The stochastic force g{t) depends on the initial 
conditions of the bath operators and is responsible for the diffusive behavior. The 
explicit but cumbersome expressions for g{t — s) and g{t) can be easily derived in 
terms of the 10 parameters 6|, ll^. Equation (|3]) can be taken as the starting point 
for Markovian and non-Markovian descriptions. The Markovian description would 
arise when there is little memory, i.e. p{t) is close to a Dirac delta and g{t) is 
close to Brownian motion. Non-Markovian dynamics arise of course for kernel not 
approximative delta-like. 


We now focus on the alternative mode of description - a chain model for the 
bath - which turns out to be more suitable to study the short time behavior of the 
dynamics. A system interacting with a chain of hrst-neighbor interacting harmonic 
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oscillators is described by the following Hamiltonian: 


H, 


CHAIN 


P 

2M 


N 


N 


V{x) + DxXi + Y, Dt-iXt-uXk + Yn (Pt + Xllxl) . (4) 


k=2 


k=l 


Here x is the position of the system, a tracer particle, and are positive cou¬ 
pling constants, are the position operators of the chain oscillators, Pk are their 
conjugated momenta, and their frequencies. The Heisenberg equations for Xk 
read 

|jx(() =-r. X((), (5) 

where = (Xi,... X^), and T is the following tridiagonal matrix; 


^ Q ^ 

-Di -D2 ... 

0 -D 2 ... 

; ; ; ■.) 


( 6 ) 


We may think of dlj) as dehning a microscopic model in its own rights, but in 
this paper we wish to take as fundamental the 10 model and aim at achieving a 
description in terms of (jll). We shall refer to the chain description as if it were 
a model and call the X^ the environment. Accordingly, we want to be able to 
relate our subsequent analysis of the chain model to the 10 model at any time. In 
order to do so, we require the Hamiltonians ifjo, ifcHAm to be “equivalent”, in the 
sense that they give the same dynamics for the tracer particle x. Accordingly, the 
parameters entering i/cHAm are not free, but they have to be particular functions of 
the parameters of Hjo. Moreover, the chain oscillators X have to be specihc linear 
combinations of the independent oscillators q. Let us introduce an orthogonal N xN 
matrix O, and dehne 

Xj = '^Ojkqk. (7) 

k 

In order to satisfy our requirements, the matrix O has to be such that the equa¬ 
tions ([2]) and ([S]) are equivalent. Substituting Eq. ([7]) in ([S]) one easily Ends that 
this requirement reduces to that of determining the orthogonal matrix O such that: 


T = O . 


( 8 ) 
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The problem of determining a matrix starting from its eigenvalues is known under 
the name of Inverse Eigenvalue Problem (lEP). Kindred relations between the two 


environments have already been considered in the literature. In 


19| the authors 


define “collective modes” through a transformation that makes use of a hierarchical 
baths construction. The parameters of the chain are not obtained analytically, but 
fitted with some experimental data. A similar hierarchical transformation for an 
infinite chain is used in 20l], where the authors derive (formally) the chain parameters 


from the propagator of the initial GLE. Another approach for infinite chains has been 


proposed in 


2l| . where the collective modes are obtained exploiting the properties 


of orthogonal polynomials. 


A. Inverse Eigenvalue Problem 


lEPs arise in different areas of theoretical and applied sciences , like e.g. vibration 
theory, control theory, particle physics, geophysics, and engineering (iTI. Il8| . In 
general, an lEP consists in finding the entries of a matrix T, starting from its 
eigenvalues and some additional initial data. In our case the matrix T is a Jacobi 
matrix, i.e. positive semi-definite, symmetric, tridiagonal matrix. A N x N Jacobi 
matrix has 2N — 1 free entries, while the 10 model gives us N conditions (the 
eigenvalues Uk). In order to find a unique solution to our lEP we need N — 1 
supplementary conditions, which are obtained by exploiting the equivalence between 
Hio and TTchain- In particular, matching the interaction terms of Eqs. o and (0]) 
one finds that DXi = which by means of Eq. ([7]) implies Oik = D~^Ck- 

The knowledge of the first line of O provides enough conditions to uniquely solve the 
Jacobi lEP. The algorithm to determine O and T is standard so we do not repeat 
it here. We do refer only to the results essential to this paper, further details can 
be found in {^, 181. The entries of the matrix O read 


Ojk 


G-i \ 

\i=i ) 


(9) 


where Pj{X) is the characteristic polynomial of the j-th leading principal minor of 
T, evaluated in A. Note that the explicit expressions for the Pj are determined 
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recursively exploiting the following recurrence relation: 

Pj+i(A) = (Jif - A)Pj(A) - (10) 

with P_i = 0. Once the transformation matrix O is determined, the entries of T 
are given by the following relations: 


k 

(11) 

Dj = — ^ ^ ujj^OjkOj^ik ■ 

(12) 


k 


These equations complete the set of parameters that we will need in the following 
discussion. From now on we will consider the matrix T as known, i.e. as fully 
determined in terms of the parameters of the 10 model. The dynamics given by 
Eq. (j3]) with the parameters here dehned is equivalent to that given by Eq. o. 
Moreover, note that the VLj are not necessarily ordered with respect to the index 
j. The chain oscillators X manifestly represent a particular collective behavior of 
the independent oscillators. Unfortunately, their structure in terms of q is rather 
cumbersome and does not allow for a straightforward physical interpretation. 


III. NON-MARKOVIAN DYNAMICS OF AN HARMONIC OSCILLATOR 

The aim of this section is to determine how the collective modes X affect the 
dynamics of the system. For this we need to solve the set of equations of motion (|5]). 
In order to do so, we approximate the tracer particle potential V{x) (cf. (0])) har¬ 
monically, i.e. we consider now an harmonic oscillator with proper frequency hi. The 
Heisenberg equations of the open system can then be explicitly written as follows: 

= -VP‘x{t) + DXi{t) (13) 

= -n^Xiit) + + AA+i(t), 1 < * < A - 1 (14) 


(15) 
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We rewrite these equations in integral form; 

At) =h+ f (16) 

Jo “0 

~ (A-iX,_i(s) +AX,+i(g))cfg , 1 < * < iV-l(17) 

Jo 

=/at + [ ^ ^ - ^—DN-iXN-i{s)ds, (18) 

Jo ^ W 

where we have relabeled hlo = = x, and 

Mt) = W(0) cos[Qd] + . (19) 


In order to obtain the dynamics of x{t) in terms of the Xj, we substitute recursively 
Eq. (ITT)) in Eq. (fTUj) . One can prove that substituting the equations up to X„(f) (for 
every n < N), x(t) reads 


n / 2 


it) = fn{t) + XI I n ^ “ s)Xi_i{s)ds 

i=l \l=o “v * 


n 

1=0 


D, ' 


Q 


I / Jo 




( 20 ) 


where recursively 


Ki{t — s) = / i^j_i(f —/) sin[r2j(Z — s)](iZ 


( 21 ) 


(ng)/ 


- s)fi{s)ds, 


with Ko{t—l) = sin[r2(f—/)], /o(t) = fo{t), and /i(s) given by Eq. (IT^ . Equation fOUj) 
is obtained by induction. To see its correctness for n = 1 substitute Eq. (ITT)) for 
Xi(f) in Eq. (ITbl) to hnd 

x{t) = fo{t) + ^ ^ sin[0(f - s)]fi{s)ds 

rt DD 

mT / ^i(^-5)-^2(s)ds. (22) 

1 «/ 0 1 «y 0 


This equation can be easily recast in form of Eq. (|2n|) . Assume now that Eq. (|2ni) 
is true for a generic n < N, and substitute Eq. (1T71) for X„+i(s) in the second line 
of Eq. (120|) . After some simple manipulation one can show that Eq. (|20|) results for 
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n + 1, i.e. 


n+l 


x{t)=fn+i{t)+Y^ n 


A 1 A-i 


i=l \Z=o 


/ A 


Ki{t - s)Xi_i{s)ds 


I n ^ ) / “ s)X„+2(s)c?s 

\z=o 


(23) 

It is important to note that equation 0201) is exact. Although only n equations 
have been substituted into x{t), the second line of Eq. (|20l) contains all the informa¬ 
tion regarding the evolution of the remaining N — n modes, encoded in X„+i(s). It 
is easy to show that, if one substitutes all the N equations for X*, x{t) satisfies the 
following equation; 


N 




i=l \l=o 




Kiit — s)Xi_i{s)ds. 


(24) 


Indeed, Dzv = 0 by definition. Accordingly, the second line of Eq. (|2ni) is null. 

This equation determines the dynamics of x{t) in terms of the full set of Xj. One 
can identify two different contributions to the evolution of x\ one is given by /Ar(f) 
that collects the initial conditions of the collective modes; the second is a purely 
non-Markovian integral contribution which involves the whole past evolution of the 
collective modes. 

In order to study the evolution of x{t), recall that Xq = x and rewrite Eq. fl2Tl) 
as follows: 


x{t) = 




t) 0 

where the function Fn {t) collects all terms of Eq. 


Ki{t — s)x{s)ds + F]s[{t), (25) 

that do not depend on x: 


N 




i=2 \Z=0 


Di 


Ki(t - s)Xi_i{s)ds , 


(26) 


Equation (|25l) explicitly shows that the dynamics of x is ruled by an integral equa¬ 
tion. Furthermore, we stress that Eq. 0261) can be rewritten in terms of the indepen¬ 
dent oscillators by exploiting Eq. ( 0 : 


xA>) = E Ojk I gfc(O) coscufcs -F 4(0) 


sm ujkS 


(27) 


Recalling the definition of Eq. 
the kernel Ki [t — s) is 

K,{t-s) = 


, and performing the integration, one finds that 


rii sin[r2(f — s)] — n sin[r2i(f — s)] 


(28) 
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Since this is a simple combination of two sine functions, Eq. 
standard techniques 


can be solved using 


22|. The solution reads 


x{t) = F^it) + 




hlh2(h2 - h?) io 


(/i 2 sin[/ii(t - s)] - sin[/i 2 (f - s)]) FN{s)ds , 

(29) 


where 


hi,2 — 


112 +^ ^ A = {n^ -nj 


2\2 


4T)2 


(30) 


Note that in order to have real values for /ri ^2 and avoid multi-valued x{t), one needs 
the condition A > 0 to hold true. Equation (|2^ is the main result of this paper; it 
displays the exact solution of our problem, representing the dynamics of an harmonic 
oscillator under the influence of a chain of N harmonic oscillators. We stress that 
all the parameters entering this equation are known analytically. Furthermore, by 
means of Eqs. fl26l) - (1?ri) , one can rewrite Eq. (|2^ in terms of the qk- this implies 
that for any sets of parameters and initial conditions of the 10 model, Eq. (129|) also 
provides the exact dynamics for this model. In this sense, Eq. (|29l) describes how 
the non-Markovian dynamics of the system depends on the microscopic motion of 
the 10 bath constituents. 

The dynamics of x{t) is determined by the function Fn {t) dehned in Eq. 0261) . 
F 7 v(t) is stochastic since it depends on the initial conditions qk{0) of the environmen¬ 
tal oscillators which are unknown. We then see that, as expected, x{t) displays the 
same diffusive behavior as the one given by stochastic force g{t) in Eq. ([3]). More¬ 
over, Fi^{t) displays two different contributions: the one given by /Ar(t) describes 
the behavior of the effective modes as if they where free oscillators. The second 
contribution depends on the interaction among the modes and the dynamics of the 
full chain. An important feature of Eq. (|2^ is the way the peculiar environmental 
structure and the non-Markovian features show up; on the one side each A* con¬ 
tributes to the dynamics via the integral kernel Aj+i, reflecting the physical ordering 
among the effective modes. On the other side, the kernels themselves have a nested 
structure [see Eq. fl?Il) ]. It is indeed this particular structure that paves the way 
for obtaining an effective description of the non-Markovian dynamics. One should 
also mention that such a structure is similar to the one obtained by Mori with his 


continued fraction description of the generalized Brownian motion 


10 


, 123]. This 


suggests that our derivation is a particular (and exactly solvable) case of the general 
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projective technique proposed by Mori. However, we could not provide any direct 
evidence for this. 

As previously mentioned, the GLE © is the reference result for non-Markovian 
dynamics. It is then essential to compare our result to the GLE and show which are 
the advantages given by the microscopic description. One major difference is that 
the GLE is an integro-differential equation that needs to be solved, and remark¬ 
ably the solution method (and its succesfullness) strongly depends on the spectral 
properties of the bath. In other words, the solution cannot be easily obtained for 
every set of environmental parameters. On the contrary, Eq. fl2U|) is the solution of 
Eq. (1231) . and it explicitly displays the dynamics of the system for any environment. 
Furthermore, since by construction the dynamics of the 10 and chain models are 
equivalent, Eq. (|2^ is also the solution of the GLE ([3]) for an harmonic oscillator. 
This result gives an important contribution to the held as it allows for the exact 
treatment of a wider class of environments. Moreover, besides extending the range 
of applicability of the GLE, Eq. has also the value of being easier to analyze both 
at the analytical and numerical levels. It is worth mentioning that while in Eq. (jHj) 
one can quite easily perform the Markov limit, obtaining the Langevin equation in 
the classical regime, one cannot do this as easily in Eq. fl2^ . The GLE and the chain 
equation hence are two complementary descriptions, from which different effective 
descriptions can arise: the hrst is suitable to understand how the environment phe¬ 
nomenologically affects the system, the second to understand how non-Markovian 
features emerge from the underlying motion on short timescales. How such an ef¬ 
fective description of non-Markovian dynamics can be realized will be analyzed in 
detail in a subsequent paper. 


IV. CONCLUSIONS 

In order to provide a microscopic description of non-Markovian dynamics, we 
transformed the 10 environment in a chain fashion, by solving an inverse eigenvalue 
problem. We derived the exact equation of motion for an harmonic oscillator, in 
terms of the microscopic motion of the environmental modes. This result repre¬ 
sents a fundamental step forward in the understanding of non-Markovian dynamics. 
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Equation shows how the microscopic behavior of the bath influences the system; 
the environmental oscillators act in a specific collective way, dehning some collective 
modes. An interesting feature of the dynamics is that each collective mode Xj acts 
through a kernel Xj+i; this peculiar structure will be subject of further studies. 
Remarkably, Eq. (j2^ provides the solution of the most important phenomenological 
equation describing general open quantum systems. This result allows to shed new 
light on the understanding of non-Markovian phenomena, allowing for an easier and 
deeper analysis both at the numerical and analytical levels. 
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